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Abstract 

We consider a dark energy fluid with arbitrary sound speed and equation of state and discuss 
the effect of its clustering on the cold dark matter distribution at the non-linear level. We write 
the continuity, Euler and Poisson equations for the system in the Newtonian approximation. 
Then, using the time renormalization group method to resum perturbative corrections at all 
orders, we compute the total clustering power spectrum and matter power spectrum. At the 
linear level, a sound speed of dark energy different from that of light modifies the power spectrum 
on observationally interesting scales, such as those relevant for baryonic acoustic oscillations. We 
show that the effect of varying the sound speed of dark energy on the non-linear corrections to 
the matter power spectrum is below the per cent level, and therefore these corrections can be 
well modelled by their counterpart in cosmological scenarios with smooth dark energy. We also 
show that the non-linear effects on the matter growth index can be as large as 10-15 per cent 
for small scales. 
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1 Introduction 



The accelerated expansion of the universe, first detected with Type la supernovae observations [Tj|2] , 
is one of the major problems of cosmology. Today, there is ample evidence for it from the combination 
of supernovae data with the cosmic microwave background and baryon acoustic oscillations [Hid] ; 
and even from the cosmic microwave background alone [5] . The task of discriminating between 
different theoretical possibilities able to describe the acceleration is of the utmost importance not 
only for cosmology but also for the field of high energy physics. In this context, the ultimate question 
that we would like to answer is what is the dark energy behind the acceleration. 

It is well known that any smooth expansion history of the universe compatible with the obser- 
vations can be described in the framework of general relativity by constructing an effective time 
dependent equation of state which is given by the time evolution of the Hubble paramete 10. This 
means that a theory of modified gravity cannot be told apart from a model based on general rel- 
ativity on the single basis of the determination of the background dynamics, even if done at very 
high precision. 

The expansion history also affects the way in which cold dark matter inhomogeneities grow 
to form structures in the universe. Clearly, a fast expansion slows the growth rate. The simplest 
equation to approach to this phenomenon in general relativity can be written, using conformal time, 
as: 



The effect of the background on the matter perturbation S m is due to the time dependencies of the 
conformal Hubble parameter "H and the relative matter density £l m . Therefore, the evolution of 
inhomogeneities provides an alternative tool to gain insight about the properties of dark energy. 

Since we know next to nothing about the actual physics causing the acceleration of the universe, 
we can adopt a phenomenological point of view and effectively describe it as a fluid in general 
relativity, as we do for matter. Such a fluid will also develop perturbations, as matter does. The 
time evolution of these perturbations will depend on the properties of the fluid: its equation of 
state, sound speed and anisotropic stress. The growth described by (jl.ip is then only valid in the 
case in which the dark energy fluid represents a cosmological constant. For any other case, the 
equation needs to be modified to take into account the effect of dark energy perturbations. These 
fluctuations change the right hand side of (jl.ip , modify the metric perturbations as well and obey 
another second order dynamical equation which also depends on the matter ones [7] • It is unclear 
whether observations related to the perturbations can tell us if modified gravity or scalar field 
models within general relativity are to be preferred, but they are undoubtedly a great source of 
information to learn about the properties of the acceleration. 

Most theoretical studies related to dark energy fluctuations that have been done so far are in 
the framework of linear perturbation theory [71424] and so are all the works on present constraints 
and forecasts for detection , but given the increasing precision of present and future planned 

probes, the relevance of non-linear effects has to be taken into account. In particular, the power 
spectrum of matter fluctuations is sensibly affected by non-linearities for wavenumbers larger than 
0.1/iMpc -1 , where a good deal of information can be gained from galaxy surveys. In this work, 
we asses the importance of these effects for the growth of matter fluctuations in the presence of a 

^See for instance [6] for an explicit expression. 
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dark energy fluid. We focus on the effect of the sound speed, assuming that the dark energy has no 
anisotropic stress. 

Presently, there are no observational constraints on the sound speed of dark energy because the 
available data are not powerful enough. However, it has been shown that combining the cosmic 
microwave background with large scale structure data it will be possible to set a lower bound on it 
in a near future [27,37,38,40] . The precise bounds will be affected by the effect of the non-linear 
evolution of the perturbations. 

In [3T] a "quasi-nonlinear" approach was considered to study the spherical collapse approxi- 
mation with dark energy fluctuations with a small sound speed. Maximal clustering of the dark 
energy perturbations occurs if their sound speed is zero. The spherical collapse was analysed for 
such a case in [42] . Second order solutions for the energy density and velocity perturbations in 
Eulerian perturbation theory were found, also for a vanishing sound speed, in |43] . The lowest 
order (tree-level) contributions to the bispectrum were also calculated in the same work. 

Dark energy fluctuations are typically suppressed with respect to those of matter and this makes 
their detection difficult. The change that they produce on the matter growth S m /a with respect to 
what is expected in a smooth wCDM cosmology without dark energy perturbations can be estimated 
from the growth index 7 [44] , which is modified by ~ 5% at most [7] . The matter growth itself, 
normalized at z —> 00 and evaluated today, changes by at most 1% , depending on the scale. These 
small changes are however large enough to be comparable to the predicted sensitivities and hence 
one could hope to connect with dark energy perturbations by measuring the matter growth. These 
results are based on linear perturbation theory and therefore it makes sense to see how they can be 
modified when non-linear corrections are taken into account. We will see that indeed the non-linear 
effects greatly affect (to the level of 10-15 %) the matter growth index for small scales and therefore 
they should be accounted for when the index is used (for model comparison, for instance). 

We concentrate on scales below the Hubble horizon, concretely on the relevant ones for the 
observation of the matter power spectrum with large scale structure surveys, and use the time 
renormalization group (TRG) [35] , which is a semi-analytical method to solve the Euler, continuity 
and Poisson equations resumming classes of perturbative corrections at all orders. In this work, 
we set a useful theoretical framework to deal with non-linearities in the Newtonian approximation 
for a mixture of matter and a non-pressureless fluid characterized by a non-adiabatic sound speed. 
We write the non-linear continuity and Euler equations for subhorizon scales and build a consistent 
approximation of the non-linear terms that allows the application of the TRG. We also explain the 
relevance of the sound speed and where it comes into play. While in the case of pure pressureless 
quintessence, the sound horizon of dark energy has zero length and fluctuations cluster for all scales 
smaller than ; if the sound speed is different from zero, scale dependent effects become relevant 
already at the linear level. When the sound speed is non-zero, at this level, the time variation of 
the pressure of dark energy becomes important and, in consequence, the velocity perturbations of 
dark matter and dark energy are not equal (the two fluids are not comoving), in contrast to the 
case of pressureless quintessence analysed in [43J . 

We apply the TRG to obtain the non-linear power spectrum of density fluctuations. Since 
at the linear level dark energy perturbations are suppressed with respect to the matter ones by 
a 1 + w factor, a non-linear treatment is strictly needed only for matter perturbations. This 
observation simplifies the implementation of the TRG approach considerably and is at the basis 
of the approximation that will be described in Section [6]. Our results indicate that the hierarchy 
between the two components that occurs at the linear level is maintained at the non-linear one, 



2 



confirming the consistency of the approximation. 

As already pointed out in |42t ,43j , if the dark energy component clusters on observable scales, 
it may become problematic to disentangle its effects from those of cold dark matter, e.g. in lensing 
measurements, so that a case by case discussion of the different observables is needed. Therefore, 
it makes sense to consider also the power spectrum of the full clustering fluid, i.e. cold dark matter 
plus dark energy. We comment on this issue examining the growth function for the matter and full 
spectra at the linear and non-linear levels. 

The structure of the paper is as follows. In section [2] we describe the Newtonian approximation 
for a fluid with pressure, writing the Euler, continuity and Poisson equations and checking their 
consistency with linear perturbation theory in general relativity. We show that no pressure terms 
appear in the Poisson equation and explain the importance of time derivatives that have been 
neglected in earlier works. In Section [3] we present the non-linear perturbation equations, which 
we provide to the level that we will need for our numerical computations. Section [4] is devoted 
to explain in more depth the general approximations that we have adopted. We review the TRG 
method and particularize it for our case of interest in Section [5l The next section, [6], describes 
the approximation that we use to reduce the number of equations in the TRG for our numerical 
solution, which is given in Section [7J. Section [8] corresponds to the conclusions. The Appendix lAl 
discusses the formulas of j43] in the context of the equations of Section [2j In our work we also use 
these formulas in Section [7J for the purpose of comparison. 

2 The Newtonian approximation 

We consider a system composed of two cosmological fluids. One of them is pressureless and repre- 
sents baryons and dark matter together. The other fluid corresponds to the dark energy component 
of the universe. Since we are interested in studying the non-linear effects of dark energy on the cold 
dark matter power spectrum and we focus on the epochs of matter and dark energy domination, 
we do not need to include other components such as radiation. We will label quantities relative to 
the matter component with a subindex and those referring to dark energy will be tagged with 
an V . Whenever that an equation involving fluid quantities without labels appears in the text, it 
will apply to any of the two components. 

In this section we will present the equations that describe the system in the Newtonian ap- 
proximation. We will specify the properties of the fluids and describe our approximations and their 
range of validity, which will be further explained in Section [4]. Although in the following we keep the 
non-linearities in the dark energy perturbations, in our numerical calculations, described in Section 
[6], we will treat dark energy perturbations linearly, an approximation motivated by equation (|2,32p 
below. 

2.1 Cosmological fluid dynamics with pressure 

Our starting point is the following couple of equations: 





u + Hu + (u • V) u + 



VP + uP 

P + P 



+ V<^v = 



(2.2) 
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The density and pressure of the fluid are denoted by p and P respectively and the peculiar velocity 
is represented by u . Overdots denote derivatives with respect to conformal time, T~L is the conformal 
Hubble factor and V is the three-dimensional gradient for comoving spatial coordinates. These two 
equations are meant to improve the classical continuity and Euler equations [36] in their application 
to describe a fluid with pressure in an expanding cosmological background. The continuity equation 
(|2.ip differs from the classical version of |46] because it includes the pressure inside the divergence 
term and in the time derivative of the first term. Besides, the Euler equation (12. 2p also contains 
the time derivative of the pressure, which was not originally present either [36] . The quantity 4>N 
represents the Newtonian potential and encodes the effect of the metric perturbations. 

The continuity and Euler equations describe local energy and momentum conservation. They 
can be obtained from general relativity applying the covariant conservation equations |47j 

T^. M = (2.3) 

to an energy-momentum tensor of the form 

T^ u = {p + P) u»u v + Pg» v , (2.4) 

assuming that metric perturbations are small (and hence c^jv <C 1) and that the peculiar velocity is 
much smaller than the speed of light: |u| <C 1 . Concretely, these equations can be obtained writing 
the four-velocity of the fluid as = (1, u) in (|2.4p (since the Lorentz factor 7 ~ 1) and neglecting 
high order metric corrections. These lines were followed in |43j neglecting the time derivative of 
u 2 P in the continuity equation^ . However, in order to do a fully consistent treatment of the 
perturbations, this term has to be included. 

If we expand this derivative using (|2.2p to replace u in (|2.ip and we keep the terms that contain 
up to two perturbations, we obtain: 

p + 3H(p + P) + V • (pu) + P V • u + ^— ^ u • (VP + uP) - 2P u • (Tin + V(p N ) = (2.5) 

p + P 

This equation shows explicitly that there are terms of second order in the perturbations that cannot 
be neglected by throwing away the u 2 P term in (|2.ip . Notice also that, even if |u| <C 1 , it is not 
possible to neglect the term uP in (|2.2p with the argument that for small velocities VP 3> uP . 
This is because uP/(p + P) contains (after applying the continuity equation for the background) a 
term — 3wT~Lu that is of the same order of magnitude as Hu , which is the second term of ()2.2p . 

If we had also included the effect of the Lorentz factor 7 we would have got extra non-linear 
corrections proportional to u • u , which we neglect. These corrections can be neglected in our 
treatment because dark energy perturbations are suppressed with respect to those of matter (as we 
will later explain in more detail). Otherwise, these corrections should be taken into consideration. 

We will split the total density and pressure according to 

p = p(l + 5) (2.6) 
P = P(l + n), (2.7) 

where the quantities with an overbar refer to the background (time dependent) homogeneous cos- 
mology and we define the perturbations 

5p = p5 (2.8) 
<5P = Pn, (2.9) 



^This term was also absent in [42] and Hi] 
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which in general are functions of time and space coordinates. The Newtonian potential <pN can be 
expressed in terms of the total density and velocity perturbations through the Poisson equation 

V 2 <p N = 4TrGa 2 Y, {Spa-3(p a + P a )H Xa ) (2.10) 

a = m, x 

Here, G is Newton's gravitational constant, a represents the scale factor of the universe, the index 
a sums over the contributions of matter and dark energy and we denote by \ the velocity potential 
of a fluid: 

u = V X (2.11) 

Therefore, we are assuming that both components, matter and dark energy, are irrotational. This 
assumption is well justified for matter because the curl of its peculiar velocity decays as a -1 at the 
linear level. For a fluid with pressure, if we define w = V x u , we have w = (— 1 + 3w)7iw, which 
implies that the decay is even faster for negative w, also at the linear level. Then, for simplicity, we 
will assume that we have zero dark energy vorticity at all times. 

It is worth noticing that the Poisson equation (I2.10P does not contain any pressure terms, in 
contrast with what was written before in other related works |41H43j and also much earlier in [J6] . 
What remains of this section is devoted to check the consistency of the equations with linear 
perturbation theory in general relativity, using the synchronous and conformal Newtonian gauges 
for illustration. In Subsection 12.31 we will review the concept of non-adiabatic sound speed and 
emphasize the need for using its rest frame value. 



2.2 The Poisson equation 

Let us start with the Poisson equation (|2.10j) . In the conformal Newtonian gauge, within the context 
of linear perturbation theory in general relativity, the Einstein equations in Fourier space (see |48j . 
for instance) imply that the two metric scalar degrees of freedom are equal to each other and given 
by: 

k 2 4> = -4nGa 2 £ Up® + 3 (p a + P a ) (2.12) 

for any system of fluids with zero anisotropic stress. Here, the labels refer to quantities in the 
conformal Newtonian gauge, k is the wavenumber and 9 denotes the velocity divergence, defined as 
(p + P)0 = ik^5T°j , and which we can link to u in the Newtonian approximation by 

9 = V-u (2.13) 

It is straightforward to obtain fl2. 121) from the perturbed Einstein equations 

5G„ V = 8itG5T, iu (2.14) 

Following the notation of [38] , we write the FRW metric for scalar perturbations in the Conformal 
Newtonian gauge as 

ds 2 = a 2 (r) {-(1 + 2ip)dT 2 + (1 - 2</>)dx i dx i } (2.15) 
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If we express the traceless component of T l j as £ i j = T l j — |<5 i jT\ and define the anisotropic stress 
(p + P)a = —(kikj — ^5ij)T> z j , where the hats stand for unit length vectors, the Einstein equations 
read 148 : 




where the terms on the right hand side have to be understood as the sums for all the fluid components 
(dark matter and dark energy in our case). Combining (|2.16p and (|2.17[) one obtains precisely (I2.12p . 
Then, under the assumption that a = (no anisotropic stress) we have 4> = ip . Therefore, in order 
to have a pressure component in the Poisson equation, the fluid would have to have non-zero 
anisotropic stress (which is a gauge invariant quantity). Notice also that during matter domination 
T~L = 2/r and therefore 21~i + H 2 is zero. This clearly indicates that the pressure perturbation 
in (|2.18p is of the order of the terms involving the time derivatives of the potential, which in the 
Newtonian approximation can be neglected. 

Using the Friedmann equation for T~L we can express (|2.12p as follows: 



The equation ()2.10|) is the counterpart of (|2.20p in real space. We will later argue that the only 
relevant term of the Poisson equation in the Newtonian approximation is the one that depends 
on the density perturbations 5 a . But before that, we are going to introduce the parameters that 
characterize the properties of the dark energy fluid. 

2.3 Equation of state and sound speed 

We will now introduce the parameters that relate the pressure and the energy density in the fluids 
we consider. We assume that the background pressure and density of each fluid are linked through 



These two expressions define the equation of state parameter w and the sound speed c s for each 
fluid. By construction, the equation of state depends only on time, while the speed of sound may 
in principle also vary in space. If we imposed that the sound speed is solely a function of time, the 
fluid would be instantaneously barotropic at the background and perturbation levels. Notice that 
the relative density and pressure perturbations are related through 




(2.20) 



P = wp 



(2.21) 



and, similarly, we write the following relation for their perturbations 



5P = c 2 5p 



(2.22) 



c s 2 5 = wll 



(2.23) 
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Our matter fluid will be assumed to have exactly zero w and sound speed, while for dark energy 
we will allow these quantities to take constant values and assume w < — 1/3 . For our numerical 
treatment of the non-linear effects, the assumption that w and c 2 are constant is the simplest 
choice one can make in absence of much information about the nature of dark energy and it is 
a phenomenologically reasonable one since in many field models both quantities are very slowly 
varying during the last stages of the evolution of the universe. Lacking a specific preferred model 
of dark energy, working with constant parameters is a simple and justified choice. 

We restrict our analysis to fluids that interact only gravitationally and therefore satisfy a con- 
servation equation of the form 

p + 3(l + w)Up = 0, (2.24) 
which is nothing else than the background (linear) part of f)2. 1 j) . In consequence, 

w = 3(1 + w)H (w - c 2 ) , (2.25) 

where we define the adiabatic sound speed of the fluid to be given by 

P = c a 2 p (2.26) 

Clearly, c 2 = w if the equation of state parameter is constant. In the general case that c s ^ c a , it 
is said that the fluid is not adiabatic. For an interpretation of the non-adiabatic sound speed c s in 
terms of entropy perturbations we refer the reader to [25] . 

In the framework of linear perturbation theory in general relativity, the sound speed defined in 
(I2.22j) is a gauge dependent quantity. Since we are interested in studying its effect in the clustering 
of matter, it is convenient to replace c s by a gauge independent quantity that can be potentially 
measured. This can be achieved by using the following relation in Fourier space |25j : 

(c s 2 -c 2 s )5 = 3(1 + w) {c 2 - c s 2 ) p0 , (2.27) 

which (in this form) is valid in the conformal Newtonian and synchronous gauges at linear level 
and gives us the rest frame (gauge independent) sound speed c 2 . In order to do a fully consistent 
second order treatment of dark energy perturbations we would have to modify the equation (|2.27[) 
to include second order terms. However, since we will treat these perturbations linearly in our 
numerical analysis (and all the non-linearities will be in the matter component), the equation 
(I2.27P will be enough for our purposes. 

Let us finally remark that in the case of negative w , the sound speed of the fluid has to be non- 
adiabatic in order to avoid pathological instabilities. For a detailed study on this issue see |16| . In 
this work we will assume < c 2 < 1 . As mentioned in the Introduction, using cosmic microwave 
background and large scale structure data, it will be possible to set a lower bound on the sound 
speed of dark energy in a near future [271[371[38lllQ] ■ 

2.4 Linear perturbations and consistency with general relativity 

As we have mentioned earlier, at the level of linear perturbations, the results obtained from the 
Newtonian approximation should match those of general relativity for scales smaller than ~HT X . 
Before solving our full Newtonian equations (|2.5p and (|2.2p , we should check that this is indeed the 
case. To do the matching between the Newtonian approximation and general relativity we work in 
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Fourier space and use c s 2 inside the Euler and continuity equations. Expressing the energy density 
and the pressure in (|2.ip by (|2.6|) and (|2.7p , and then making use of (|2.27p to introduce the rest 
frame sound speed, we get that the linear parts of these equations for the case of constant w are 

S=-(l + w )0 - 3 (c s 2 -w)H5- 9(1 + w) (c s 2 - w) -^9 (2.28) 

c 2 k 2 

= -h- 34 2 ) <ue + — 6 + k 2 <t> N (2.29) 

l + w 

These equations look very similar to the ones of general relativity for linear perturbations in the 
conformal Newtonian gauge: 

5^ = -(i + w ) (6® - 30J - 3 (c s 2 - w) H6® - 9(1 + w) (c s 2 - w) ^9^ (2.30) 

0(c) = -h- 3c s 2 ) H9® + ^-8® + k 2 ^ (2.31) 

l + w 

Already from here, we can immediately guess that in order to get the Newtonian approximation 
we need to identify cf> and (J)n and be able to neglect the contribution of (p in 112. 30 j) . Indeed, using 
(|2.20p we see that the contribution of <f> in (|2.30p gives contributions in S and 9 that are suppressed 
by powers of H 2 /k 2 and H 4 /k 4 , respectively. 

It was found in [37] that the relative energy densities of matter and dark energy in the regime 
of pure matter domination (£l m = 1) are related by 

5 x = {l + w) 1 - 2Cs 2 5 m , n s ^k^n (2.32) 
1 — 3w + c s 



for scales smaller than 7i~ 1 (where the density perturbations are approximately gauge invariantH 
and larger than the sound horizon Hj 1 = CgH" 1 . This expression can be obtained using the 
synchronous gauge by looking for the possible solutions of the second order differential equation for 
5 X and taking into account that the growing mode for cold matter perturbations evolves in time 
as 5 m cx k 2 T 2 during matter domination. The matter perturbations act as a source term for dark 
energy fluctuations and it is this source term that originates the solution (|2.32p . 

Since the Newtonian approximation applies in this range of scales, the solution (12.32|) and the 
(gauge invariant) density perturbations computed in general relativity (in any gauge) should be 
recovered from (|2.28p and (|2.29p in combination with (|2.20p . A good consistency test to check that 
this is indeed the case can be done from the second order differential equation for 5. To obtain this 
equation, first in the conformal Newtonian gauge, we differentiate the continuity equation (|2.30p 
with respect to conformal time and use the Euler equation (I2.3ip and the continuity equation itself 
in the result. In this way we can eliminate the velocity divergence 6>( c ) . What we get is an equation 
(valid at any scale) that only contains 5^ and its first and second derivatives on the left hand side; 
and the metric potential (j) and its derivatives on the right hand one: 

§(c) + [3 (£2 _ w ) n _ j-] j(c) _ 3 ^2 _„,)[(! + 3^ - 6c, 2 ) n 2 + 2UJF] 5^ + c s 2 k 2 5^ 

= (l + w) (30 - 3T<p - V(f>) , (2.33) 



f See [48] or [37], for example. 
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where 



V = kT + 9(c/ - w)W (2.34) 

T = -9 (1 + 3^) Cs V W , H 3 - (1 - 3c s 2 )^ (2.35) 

Notice that if we follow the same procedure using the equations (|2.28j) and (|2.29[) instead, we arrive 
to an expression which has formally the same functional structure as (|2.33p with the difference that 
the derivatives of 4>n are absent. This is consistent with our previous guess that the time derivatives 
of (j) must be negligible in the Newtonian limit. For scales k 3> % , the right hand side of (|2.33D 
becomes equal to 

(1 + w) (30 + 3(1- 3c s 2 ) H<j) - k 2 ^j = (1 + w) (5 m {c) + m'J c) - 9c 2 H^j , k~>U (2.36) 

If we had done analogous operations working in the synchronous gauge, we would have obtained 
the same functional form for the left hand side of (12.33P (as explained in 0E7]) and the following 
sum 

(1 + w) UJ S) + (1 - 3c s 2 )n5 m (s) ) , k^n (2.37) 



(s) 

at the right hand side, where 5 m stands for the matter density perturbation in the synchronous 
gauge. Taking into account that the gauge differences of dark energy and dark matter perturbations 
are related by 

#)-#> = (! + «,)($)-#)) (2.38) 
and expressing the time derivative of the potential as 

<P = ^ («# " + \S® (2-39) 

one can check that the second order differential equation in the synchronous gauge [37] is equal 
to (|2.33p in the limit k <C H (up to negligible terms that are suppressed by factors of T~L 2 /k 2 or 
higher powers) . Therefore, both gauges do indeed give the same result for subhorizon scales, as 
expected. In addition, notice that for subhorizon scales and in matter domination 5 m = 5 m 7i and 
therefore (|2.39j) tells us that <j) — >■ and (12.32j) is also valid in the Newtonian approximation. In 
consequence, we see that the conformal Newtonian gauge is the best suited to arrive to the equations 
in the Newtonian approximation and we have checked the consistency of these with the synchronous 
gauge at small scales. From the arguments above, we can expect that working with the following 
Poisson equation should be an excellent approximation: 

-k 2 <t> N = ^H 2 Yl n ^S a , (2.40) 

a = m, x 

because the contributions of the velocity diverg ences 0^ ^ / }~Lb GL from (|2.20p to the Euler equation is 
suppressed (see also Section U]) . 

Notice also that looking at (j2.30j) or (|2.28p one could naively expect that it would be correct 
to neglect the term in 9 that is proportional to Tl 2 /k 2 in those equations. However, such an 
approximation is not consistent for early times because it fails to reproduce the result (|2.32p , as 
one can easily check. The subtle reason behind this fact is the presence of the term proportional 
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to c s 2 k 2 8 in the Euler equation. This feature is characteristic of fluids with pressure and therefore 
does not occur for pure dust. Moreover, let us remark that it is this same term in (|2.33p that allows 
the equations to be consistent in different gauges. 

We can also express the velocity divergences in terms of the matter perturbation deep in the 
epoch of matter domination (for Q m =l) . One can check, using (12.280 , that in the Newtonian 
approximation the relation is 

/ 6c s 2 (c s 2 - w)\ 

e x = -1+ / y d-)H8 m , Hs^k^H. (2.41) 

y 1 - 3w + c/ J 

Clearly, the result obtained from (|2.30p neglecting <f> (which we can do for those scales) is the same. 
This relation, together with (|2.32|) , will serve us to define the initial conditions in our numerical 
analysis. 

Having seen that our equations are consistent with general relativity not only at the background 
level (|2.24p but also for linear perturbations around it (taking the limit from two different gauges), 
we can move on to study the non-linear dynamics. 

3 Non— linear perturbations 

To find the non-linear continuity and Euler equations for the perturbations, it is convenient to write 

VP = ^V(c s 2 5) (3.1) 
P= {pc s 2 8)' + c 2 p (3.2) 

and then use (|2.27p together with the Einstein equations for zero intrinsic curvature 

3H 2 = 8ttG a 2 ^2 p a (3.3) 

a 

m = -4vrG a 2 ^2(p a + 3P a ) (3.4) 

a 

In Fourier space, the result is: 

5(k) + 3 (c s 2 - w) HS(k) + (1 + w) (l-9 (w - c 2 ) ^) 6{k) 

+ {l + c s 2 ) j d i pd 3 c i 5 D (k-p-q)a(ci,p)e(q)5(p) + 0{2) =0 (3.5) 

9(k) + (1 - 3c 2 ) ne(k) - 7^V( k ) " k ^N 

(1 + w) 

+ (l-c s 2 ) J d 3 pd 3 q5 D (k-p-c l )P(c L ,p)e(q)e(p) + O(2) = 0, (3.6) 

where 0(2) stands for the second order terms that we have not written explicitly. In order to 
compute all the contributions to those terms we would need to obtain the completion of (|2.27p at 
second order, although some contributions can be readily obtained using just (|2.27p . Our aim in 
this work is to compute the non-linear corrections to the power spectra (of matter and matter plus 
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a(q.p) = 2 — ' /nq>p) = — — • ( 3 - 7 ) 



dark energy) so in principle, all the 0(2) terms should be included. However, given that dark energy 
perturbations are suppressed with respect to those of matter, we will not need the second order 
terms for dark energy in our numerical calculations and the expressions above will be enough. This 
is justified by the expression (|2,32p , which indicates that for the values of w and c 2 that we are 
interested in, S x is roughly an order of magnitude smaller than 6 m . The functions of momenta, a 
and f3 , that appear in the explicit integrals of (|3.5p and (|3.6[) are the usual ones for dust (see for 
instance |49j ) : 

(q + p) • q o, v (p + q) 2 (q -p) 
— — , p(q,p) = tt?-^ 

q z 2q z p z 

The equations (I3.5D and (I3.6P become the standard ones for matter taking w and c s 2 to be equal 
to zero. This is so because the second order corrections 0(2) are either proportional to the sound 
speed c 2 or the equation of state w . 

4 Distance scales and non— linear power counting 

In this section we explain in more depth our approximations and comment on the relevance that 
the sound horizon of dark energy has on them. 

We are interested in scales much smaller than the Hubble distance, i.e. k S> H , which is the 
region of validity of the Newtonian approximation. However, "H" 1 is not the only relevant scale in 
the problem; apart from it we have to consider the sound horizon of dark energy, which is defined 
as TiJ 1 = CgH^ 1 . If the rest frame sound speed is equal to the speed of light, the sound horizon 
and the Hubble scale coincide. This is what happens, for instance, for quintessence models with a 
canonical kinetic term. However, in the general situation in which H s ^ % there will be two regions 
of interest for us, characterized by wave numbers k larger o smaller than T~L S . At the linear level, 
for scales smaller than the sound horizon (k > 7i s ) , dark energy perturbations display oscillations 
that do not occur outside it. 

The fact that we will focus on k ^> H at all times, allows us to make a power counting scheme 
to guess the relevance of the different terms that appear in the continuity, Euler and Poisson 
equations. It is well known that for scales smaller than the Hubble distance , 6 ~ T~L5 at the level of 
linear perturbations. This is the classical linear growth behaviour for matter perturbations, but it 
is also true for dark energy fluctuations, as can be checked neglecting the non-linear terms. 

If the relative energy densities 5 a remain smaller than 1, the non-linear terms will generically 
give small contributions. Let us first look at the continuity equation. Going to Fourier space and 
using the Poisson equation, the value of 5/% is determined by linear contributions of order & c {i) 
and non-linear ones of order A c ( n /) where 

A c(0 ~ (5) x (0 (1) + {U 2 /k 2 )) (4.1) 

and 

A c(n0 ~ (5) x (a c(1) + 0(5 a )xO (n 2 /k 2 ) xfl + O H^j ) \ , (4.2) 

where 5 a refers to any of the fluids that are present (concretely, matter and dark energy). 

Since we are interested in scales such that k 3> H we can be tempted to neglect the terms of 
order proportional to {U 2 /k 2 ) and (W 4 /A: 4 ) , but we need to be careful because we have already 
seen in Section [23] that neglecting a term of order T~L 2 /k 2 in the linear differential system would lead 
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to a wrong result in the limit of pure matter domination. It is only through a numerical solution 
of the whole system of equations that we can be certain of whether we can leave aside this kind of 
terms or not. This is actually what we have done, reaching the conclusion that they are actually 
irrelevant for the redshifts that we are interested in. 

Something similar occurs for the Euler equation. Taking the divergence of (I2.2p we see that this 
equation has the order structure 

^ ~ OW x (1 + m) x (o (*) + O ( * )) + 0( M x (l + O (f )) (4.3) 

Collecting the lowest order terms, we see that 9/T-L 2 is given at the linear level by terms of the 
order 

A E(l) ~ O(i) x ( 1 + (£) +0 (£) +0(Wx ( 1 + (£)) (4.4) 

At this level, depending on the scale of interest being below or above the sound horizon, the term 
of order 5k 2 /H 2 may be relevant or not, becoming the dominant one in the limit of very large k. 
The same occurs at the non-linear level. 



5 The time renormalization group 

The time renormalization group (TRG) is a useful technique, introduced in [35] , to compute non- 
linear corrections to the power spectrum and higher order correlators in the framework of Eulerian 
perturbation theory. The principal idea behind it is the recursive application of the equations for 
the perturbations. For any field <p(k, r) formed by several components we define: 

(^ Q (k)^(q)) = 5 D (k + q)P a6 (k) (5.1) 
(<PaQs)(p b (q)(p c {p)) = 5 D {k + q + p)B abc (k, q, p) (5.2) 
((p a (k)ip b (q)(p c (p)ip d (r)) = 5 D (k + p + q + r) Q abcd (k , q , p , r) 

+ 5 D {k + q) S D (p + r)P ab (k)P cd {p) 
+ <5 D (k + p) (q + r)P ac (k)P 6d (q) 

+ <5z,(k + r) 5 D (q + p)P ad (k)P bc (q) , (5.3) 

where we have omitted the time dependencies to abbreviate the notation. As usual, P a ^(k) is the 
power spectrum, B abc (k, q, p) the bispectrum, and Q a bcd(k , q , p , r) , the connected part of the 
four-point function, the trispectrum. We will neglect the trispectrum in our computations; and 
therefore the four-point function will be fully given by the power spectrum. 
If we define 

¥>* = e-" (5 m , -6 m /n , S x , -9 X /H) , (5.4) 
where the superindex ' ' ' simply denotes matrix transposition and 

a 

rj = log — , (5.5) 
we can write the continuity and Euler equations for matter and dark energy in the following form 

¥>a( k > V) = -^a&(k, r])ip b (k,r]) + e TI J d 3 pd 3 c[5 D (k - p - q)7 a6c (k, -p, -q, r))<p b (p, rj) <p c (q, rj) , 

(5.6) 



12 



where primes denote derivatives with respect to rj . The linear evolution of the perturbations is given 
by the matrix £l ab (k, rf) and the non-linear contributions come from the vertices -y abc (k, p, q, rj) . 
Recursive iteration of the equation (15, 6p gives 

d v (fa^fb) = -£lac{<PcPb) ~ ftbciPaPc) + e V ^acd{fcfdfb) + e V -fbcd(<fafcfd) , (5.7) 
dr, (<PaW£c) = -QadiVdPbCfc) ~ ^bdi^a^Pd^c) ~ ^cdi^a^b^d) 

+ e^ladeifdVeCPbCPc) + e V '-fbde{ l Pafdfe i Pc) + e ' ^cdei^a^bfd^e) , (5.8) 

where we have omitted the momentum dependencies because they can be tracked with the indices. 
Using in these two equations, the definitions of the correlation functions that we have introduced 
above, we get the following equations for the power spectrum and the bispectrum: 

^P o6 (k) = 

- 0, c (k)P c& (k) - Q bc (k)P ac (k) 

+ e v J d 3 q [7acd(k, -q, q - k) B hcd (k, -q, q - k) + B acd (k, -q, q - k) -f bcd (k, -q, q - k)] 

(5.9) 



d v B abc (k, -q, q - k) = 

- rt a d(k)B dbc (k, -q, q - k) - fi 6d (-q)B a(fc (k, -q, q - k) - fi^q - k)-B aM (k, -q, q - k) 
+ 2e v [y ade (K -q, q - k)P<f 6 (q)P ec (k - q) + 7&<fe(-q, q - k , k)P dc (k - q)P ea (k) 

+ 7c(fe (q - k, k, -q)P da (k)P eb (q)] , (5.10) 

where all the functions are also time dependent and we have used the vertex symmetry (|5.13p to 
group terms coming from different correlators in (|5.8p . 

The matrix that gives the linear evolution of the perturbations in the case of non-interacting 
dark matter and dark energy is 



Qab 



/ 


1 











\- 


2 a L m 


to 


(ESD 



-1 

z ^ H ^ 2 iLm Tz 



9 H 2 





■In, 



1 — 3w + 3c s 2 





-(l + w) (l + 9(c s 2 - 
H> 



\ 



w 



nr 



J 



determined by the vertices 

labc (k, p, q) = S D (k + p + q) T abc (k, p, q) , 
which, by construction, have the permutation symmetry 

r a 6c (k, p, q) = T acb (k, q, p) 
The ones that correspond to the terms that we have written explicitly in (|3.5p and ()3.6p are: 

r m (k,p,q) = -a(p,q) , 



T222(k,p,q) = /3(p,q) , 



(5.H) 
ion is 

(5.12) 
(5.13) 

(5.14) 
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which are the usual ones for matter perturbations. 

The equations (|5.9p and (|5.10p can be simplified for cosmologies for which the matrix $l ab is 
scale independent, becoming: 



d v Pab{k) 

Iacd,bef 

(k) 



ttacPcb{k) ~ Sl bc P ac (k) + e V — [I a cd,bcd(k) + hcd,acd(k)} 
^bglacd,gef(k) ~ ^eglacd,bgf {k) ~~ ^ / gl acd,beg{k) + ■^■acd,bef {k) 



(5.15) 
(5.16) 



where we define 



Iacd,bef(k) 




I dqq dpp- % cd (k,q,p) B bef (k,q,p) + (q ^ p) 

k/2 J\q-k\ 1 L J 

poo pq 

f dqq dpp-{% cd (k,q,p)[% h (k,q,p)P ge (q)P hf (p) 



(5.17) 



A a cd,bef{k) 




k/2 J\q-k 



+ %gh{q,P,k)P gf (p)P hb {k) + j fgh (p,k,q)P gb (k)P he (q)] + (?Hp)} 



(5.18) 



with 



%bc(k,q,p) = 7afe c (k,q,p)| p= _ (k+q) 



(5.19) 



and analogously for B abc (k, q,p) . 

Clearly, Q ab is not scale independent for the case of a non-vanishing sound speed of dark energy, 
but however we can actually use (|5.15p and (|5.16p reliably, simply because the error introduced by 
doing so is small. To check this, we fixed the momenta in the £l ab matrices appearing in the second 
and third terms of (|5.16p to different values below the Hubble horizon. We found that the solutions 
are nearly independent on these choices, introducing a relative error between different assignments 
that is well below 1% . 

The TRG has already been used to study not only the ACDM model [45] , but also the effects 
of massive neutrinos [50] and to take into account high order correlation functions [51] . It was also 
modified for the case of non-Gaussian initial conditions [52] , applied to study neutrino quintessence 
[53] , the matter power spectrum for a variable equation of state [5l] , the case of coupled quintessence 
[55] and dark matter haloes [56] . Here we use it to describe the non-linearities in the clustering of 
matter in the presence of dark energy perturbations. 

6 An approximate numerical solution 

Solving the full TRG of Section [5], where the field (p has four components and new vertices need to be 
introduced, requires a large computing time, even if the equations for scale free linear propagation, 
(|5.9p and (|5.10p , are employed. One can lessen the computational cost considerably by reducing 
the number of components and, in consequence, of equations. In this section we describe the 
approximation that we choose in order to achieve it. The equation (|2.32p implies that 5 X S rn 
in matter domination, but this is also true later in the evolution of the universe. As we have 
already explained, there is roughly one order of magnitude between the matter and dark energy 
perturbations and therefore S x is always linear on interesting scales, provided that 6 m remains 
sufficiently small. Neglecting non-linear terms in the dark energy fluctuations and using the Poisson 
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equation (|2.40p , we obtain the following continuity and Euler equations for matter perturbations: 



5 m (k)+0 m (k) + J d 3 pd 3 q5 D (k-p-q)a(q,p)e m (c l )5 m (p) =0 (6.1) 

e m (k) + ne m (k) + ^H 2 n m s m (1 + + J dWqfc(k - P - q)/%, p)M<i)Mp) = o , 

(6.2) 

where 6% and 5^ are computed from the linear Newtonian theory (|2.28p and ()2.29j) and we have 
done the further approximation of replacing S m by 5^ in the denominator of the third term in the 
Euler equation. This allows us to reduce the number of components of the ip field to the two matter 
variables: 

^ = e"" (5 m , -9 m /H) (6.3) 

with 

* (k ' v) = {-Pm [1 + M) I (M91 2 + f ) ' (6 - 4) 

where we have neglected the subleading terms of the Poisson equation. We have checked numer- 
ically that at late times (i.e. for small redshifts) the effect of adding those terms is irrelevant for 
scales k > 0.01 /iMpc" 1 . Notice that the only vertices that we need in this approximation: 7112 
and 7222 5 correspond to the usual ones for matter (see [45] , for instance) . This is the kind of 
approximation that was already advocated in [45] (and later used in [50] for studying the effect of 
massive neutrinos) . 



7 Numerical Results 

Here we present the results of solving the Newtonian approximation for our system of dark matter 
and clustering dark energy applying the TRG. The main observable we are interested in is the two 
point correlation function (the power spectrum). We obtain the matter power spectrum P m = (5^) 
and the total power spectrum P^(k) defined as 

p tot (k) = nl(5 rn (k) 2 ) + 2n m n x (5 m (k)5 x (k)} + nl(s x (k) 2 ) (7.1) 

Besides, we present results for the (total and matter) growth functions and matter growth index, 
which we define in Subsection 17.21 . The total power spectrum Pt t(k) is related to the total relative 
energy density perturbation of dark matter and dark energy, that we define to be: 

$tot = Qm&m + Qx^x (7.2) 

The total clustering density perturbation is dptot = Sprn + 5p x ■ This is the density perturbation that 
appears in the Poisson equation (I2.10P and therefore the one that affects the dynamics of clustering 
and structure formation. Since we want to confront this density perturbation with the total energy 
budget, we divide it by the total background density ptot = Pm + Px , obtaining the total relative 
clustering energy density perturbation (|7.2p , which is the relevant quantity in (|2.40p . In the case 
of ACDM , dark energy does not cluster because A is perfectly homogeneous and isotropic. This 
means that the total relative clustering energy density perturbation for ACDM is £1^8^ . 
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We choose the cosmological parameters to be close to those of the current best fits, assuming a 
flat universe. In particular, we take = 0.25 , Slg/i 2 = 0.0224 , h = 0.72 , n = 0.97 and a 8 = 0.8 
(defined under the assumption of w = — 1 ) . Other reasonable choices do not change our conclusions 
in any respect. The initial time, tj = , is associated to redshift Z{ n = 100 . At this time, we set 
the initial conditions for the TRG equations according to the relation between the dark energy 
perturbations and the cold dark matter ones given by the equations (|2.32j) and (|2.4ip . As input, we 
introduce the linear cold dark matter power spectrum obtained from CAMB [57] . To run the TRG 
integration we choose the initial values for the bispectra equal to zero as we neglect all the non- 
Gaussianities generated at redshifts higher than z = 100 . Recall that within our approximation the 
trispectrum remains zero at all times. 

Let us now describe the two different power spectra that we calculate. In Subsection 17.21 we plot 
the matter and total growths and explain the non-linear effect on the matter growth index. 

7.1 The power spectra 

We compute the linear P m solving the continuity (|2.28p and Euler (|2.29p equations for matter and 
dark energy, which can be expressed in terms of the matrix (|5.1ip . In practice, we work with the 
Poisson equation (|2.40j) because, as we have already explained, using the full Poisson expression for 
cpN including the velocity perturbations does not alter the results. We also calculate the matter 
power spectrum non-linearly using the TRG, applying the equations (|5.15p and (|5.16p with the 
approximation described in the previous section (in which we treat linearly the contribution from 
dark energy perturbations) . Comparing the two results we observe and describe the non-linear 
effects on the matter power spectrum. 

We proceed in the same way for the total clustering spectrum Ptot of matter and dark energy. 
The linear version of it is given by (17. ip where 5 m and S x are both computed at the linear level. 
Our approximation to its non-linear counterpart is described in Section [6]. 

In the particular case of zero sound speed of dark energy, the equations for the velocity per- 
turbations of matter and dark energy are identical. For the sake of comparison, we also compute 
the total non-linear power spectrum making use of this fact to simplify the system of equations, as 
explained in the Appendix. 

Let us notice that, strictly speaking, in order to compute a power spectrum at second order 
we should know the density perturbations at order three. This is because if we take into account 
contributions of order 5^6^ we should also include contributions of the type 5^5^ , where the 
number in parenthesis includes the perturbation order. Clearly, these two types of contributions 
are the smallest of the ones we need to consider and therefore we can work (as we do) expanding 
the continuity and Euler equations up to order two in the perturbations. 

7.1.1 Matter power spectrum 

The Figures [Hand [3] correspond to the matter power spectra for two values of the equation of state 
of dark energy: w = —0.8 and —0.9 respectively. Each curve corresponds to a power spectrum 
either linear (dashed lines) or non-linear (continuous lines). The figures show the relative difference 
between a spectrum and its linear smooth counterpart P^ , which is defined as the matter linear 
power spectrum in the limit of c s 2 = 1 and, in practice, can be calculated from (j6.4|) setting 5^ = . 
Clearly, the only parameter that affects Pg 1 ™ is the equation of state of dark energy (through the 
background expansion). This spectrum is represented in the figures by the dashed black line, it is 
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Figure 1: Relative difference between the matter spectra and their linear smooth counterpart. The spectra 
are evaluated for w = —0.8 at redshifts 2 = (top) and z = 1 (bottom). Dashed lines correspond to linear 
spectra while solid ones are non-linear ones computed using the TRG. The black lines represent the smooth 
case while the dark green (c s 2 = 0.1), grey (0.01), magenta (0.001), light green (0.0001) and blue (0) show 
the results for various sound speeds of dark energy. 
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Figure 2: Relative difference between the total power spectrum (matter and dark energy, denned as in l|7.1|) ) 
and its linear smooth counterpart. The spectra are evaluated for iv = —0.8 at redshifts 2 = (top) and 
2 = 1 (bottom). Dashed lines correspond to linear spectra while solid ones are non-linear results computed 
using the TRG. The black lines represent the smooth case while the dark green (c s 2 = 0.1), grey (0.01), 
magenta (0.001), light green (0.0001) and blue (0) show the results for various sound speeds of dark energy. 
The black dotted line corresponds to the ACDM cosmology and the orange dash-dotted one to the numerical 
approximation of the Appendix for c s 2 = . 
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scale independent and, of course, sits at zero. The black continuous line shows how P s ™ gets modified 
by the non-linear (pure cold dark matter) corrections to the Euler and continuity equations. We 
can see that for large scales the difference is minimal but at small ones it grows very fast. For 
k = 0.1 ft,Mpc _1 it is of order 6% for redshift zero in the case of w = —0.8 . 

The different colors of the curves distinguish between values of the sound speed of dark energy. 
The differences between the dashed lines (linear spectra) in each of the Figures Q] and [3] are due to 
the enhancement of the matter power spectrum that the dark energy clustering induces at large 
scales. Clearly, the closest is c 2 to the speed of light the smaller is the effect, which is maximal 
for c s 2 = 0. For scales larger than the sound horizon of dark energy (A; < T-L s ) the dark energy 
fluctuations cluster. This contributes to enhance the gravitational potential, which in turn produces 
an increment of the matter perturbations. The gentle decays that are observed at different scales for 
different speeds of sound are due to the transition between the regimes k <% s and k > 7~L S . This 
effect is similar to that of neutrino free-streaming. It can also be observed in the non-linear spectra 
(continuous lines) for large scales, but then, once the non-linearities become important (they very 
clearly do so at around 0.07/iMpc" 1 ) , the linear suppression of the spectra is washed out and more 
difficult to see. 

The feature that is observed at 0.1/iMpc -1 in Figure Q] is due to the effect of non-linearities on 
the valley of the baryon acoustic oscillations that appears to the right of the first peak (which is 
roughly located at 0.07/iMpc -1 ) . Similarly, the small bump that appears between 0.04/iMpc" 1 and 
0.06/i Mpc" 1 in Fi gures [T] and [3] corresponds to the same effect on the valley to the left of the first 
peak. One can understand these features realizing that we plot (and so the valleys become 

peaks) modified by a factor that measures the coupling between modes of different momenta in the 
non-linear power spectrum. 

A remarkable property of both figures is that the non-linearities seem roughly independent on 
the sound speed. In reality, the correction is c s 2 dependent but it varies at most by approximately 
0.5% changing the sound speed between and 1 (for the smallest scales ~ 0.1/iMpc -1 in the case 
of z = and w = —0.8 ) . The TRG resums corrections at all orders in perturbation theory. The 
leading ones correspond to 1-loop diagrams that involve the product of two linear power spectra (as 
can be seen from the pair of equations (|5.9p and fl5. lQf) ) in such a way that the order of the resulting 
corrections can be estimated in terms of the correction to the power spectrum itself. Concretely, the 
difference between the linear matter power spectra P m and P, 1 ™ is approximately P ~ (1 + A)P S '™ , 
where A ~ 2% (at most) for c s 2 = (see Figure [1]) . Therefore, the leading non-linear corrections 
from the TRG are of the order f P 2 dk ~ (1 + 2A) f dk (P^) , where f dk formally denotes 
the momentum integrals. We see from Figure [1] that P/P, 1 ™ ~ 1.1 for k ~ 0.15 /iMpc -1 (which 
roughly corresponds to the limit of validity of the TRG method). From here we can estimate the 
maximum non-linear correction due to the sound speed (happening for c 2 = 0) which turns out 
to be approximately the 0.5% difference that we mentioned above . Clearly, for sound speeds larger 
than zero the effect is even smaller and so we see that indeed the effect of the sound speed on the 
non-linearities is very small. 

The importance of the non-linearities increases as we reduce the redshift. For instance, if we 
look at the case w = —0.8 and k = 0.1 /iMpc -1 , we see that the 2.4% difference with respect to 
P s ™ observed in the bottom panel (z = 1) of Figure CD becomes, as we already said, 6% for redshift 
(top panel) . 
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k (h/Mpc) 

Figure 3: Relative difference between the matter spectra and their linear smooth counterpart. The spectra 
are evaluated for w = —0.9 at redshifts z = (top) and z = 1 (bottom). Dashed lines correspond to linear 
spectra while solid ones are non-linear ones computed using the TRG. The black lines represent the smooth 
case while the dark green (c s 2 — 0.1), grey (0.01), magenta (0.001), light green (0.0001) and blue (0) show 
the results for various sound speeds of dark energy. 
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Figure 4: Relative difference between the total power spectrum (matter and dark energy, denned as in Q7.1jl ) 
and its linear smooth counterpart. The spectra are evaluated for w = —0.9 at redshifts z — (top) and 
z = 1 (bottom). Dashed lines correspond to linear spectra while solid ones are non-linear ones computed 
using the TRG. The black lines represent the smooth case while the dark green (c s 2 = 0.1), grey (0.01), 
magenta (0.001), light green (0.0001) and blue (0) show the results for various sound speeds of dark energy. 
The black dotted line corresponds to the ACDM cosmology and the orange dash-dotted one to the numerical 
approximation of the Appendix for c s 2 = . 
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7.1.2 Total power spectrum 



The Figures [5] and H] correspond to the total power spectrum for w = —0.8 and —0.9 respectively. 
As for the matter ones, each curve corresponds to a power spectrum either linear (dashed lines) or 
non-linear (continuous lines) and we use the same color code as before. Now, we denote by 
the linear total power spectrum for c s 2 = 1 and, in practice, it can be calculated from the linear 
evolution equations setting to zero the dark energy fluctuations. We see, according to the definition 
(|7.ip , that this spectrum has now an prefactor that does not appear in the pure matter case. 
As we already mentioned before, the total power spectrum is computed from (I7.ip as explained in 
Section [6]. 

There are two extra lines with respect to the matter plots. The orange dash-dotted corresponds 
to the approximation of the Appendix and the black dotted one to the ACDM cosmology. This last 
one is out of the figures for redshift 1 due to the enhancement from O m . 

Following the discussion in Section [6], the orange line serves us to confirm the validity of our 
approximation of treating linearly the dark energy density perturbation 5 X . The fact that the differ- 
ence between the continuous blue line (non-linear, c s 2 = ) and the orange one (the approximation 
of the Appendix), that corresponds to the single fluid approximation, is so small for any scale 
means that for sound speeds different from zero, treating 5 X linearly will be even more justified. 
The equation f)2.32[) shows that 5 X < S m (during pure matter domination, but this will also happen 
at later epochs). We see that 5 X can be even an order of magnitude smaller than 5 m because of 
the (1+w) factor, therefore the linear approximation can be used for 5 X even when 5 m must be 
treated non-linearly. Besides, the clustering (and hence 5 X ) is smaller for larger values of c s 2 and 
therefore, since the plot shows that using the linear equations for the dark energy perturbations 
works well, we conclude that the linear approximation will be even better for c s 2 closer to the speed 
of light. Although it can not be read directly from the figures, the percentage difference between 
the spectra calculated as in the Appendix and with our approximation in the case of c s 2 = only 
becomes bigger than 1% for the largest equation of state that we consider (w = —0.8) , z=0 and 
k > 0.16 /iMpc -1 , at the edge of the range of validity of the TRG. 

Notice that the feature that appears in the Figures [T] and [3] at k ~ 0.1 /iMpc" 1 is present in the 
Figures [2] and 0] but stretched. The bump at around k > 0.05 /iMpc" 1 in the matter power spectra 
can also still be (weakly) seen in Figure H] at k > 0.16 /iMpc -1 , but washed out by the effect of Q x 
and the choice of a linear horizontal axis. 

Looking at the figures [2] and |4] for redshift zero, we can see that the difference between the 
relative ACDM and zero sound speed cases depends on the scale. Besides, the redshift dependence 
is a most interesting feature because it could be a way of detecting a deviation from ACDM , using 
measurements at a fixed (or several) scales for different redshifts (between and 1, or even a larger 
range) . The idea also applies to the other sound speeds and so this could be a characteristic that 
may in principle help to discriminate between different values of it. 



7.2 Growth of perturbations 

The growth of matter perturbations 

9m = — , (7.3) 
a 

introduced in [58] , is an interesting function to test, for instance, deviations from a ACDM cosmol- 
ogy due to an equation of state different from — 1 . In the limit of pure matter domination g m — > 1 . 
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It is well known that it is possible to fit g m using a simple parameterization [33] that defines the 
matter growth index j m and depends on the relative energy density of matter ; 



A good approximation for the matter growth index in the case of wCDM cosmologies is [44] 



When dark energy perturbations are taken into account (and let us recall that one needs to do so 
if w 7^ — 1) , the growth and the growth index become scale dependent and vary with the sound 
speed of dark energy and its equation of state [7] . We can write the equation (|7.4|) in a useful way 
to study these dependencies on the growth index: 



where r] is defined in (I5.5P . 

Usually, the growth and growth index are employed at the linear level but it is clear that both 
quantities receive contributions from the non-linear terms in the evolution equations. In this work 
we compute these effects taking into account the presence of dark energy perturbations. 

In an analogous way to the one that is customarily done for matter, we can define the total 
growth function gtot that corresponds to the full clustering perturbation of dark matter and dark 
energy (|7.2[) . In order to do this, we simply replace 5 m in (|7.3|) by Stat = Qm<>m + Qx$x ■ This total 
growth is relevant for structure formation in situations where the dark energy component clusters. 
The linear total growth index has been studied in [33] for the case of zero sound speed. 

Let us now focus on the cold dark matter growth function which is shown in Figure [5]. The 
smooth linear case is represented by the black dashed curve, which is indistinguishable in the graph 
from the dotted magenta (c s 2 = 0.001 linear). Let us recall that in the limit c s 2 -> 1 we recover 
the smooth dark energy case. The right panel of the Figure [5] shows that the effect of varying the 
sound speed of dark energy is mainly a linear effect, because it is essentially the same for the linear 
and non-linear curves. This effect is indeed due to the c s 2 k 2 term of the (linear) Euler equation. 
There we see that the blue line (either linear or not) is always above the corresponding magenta 
one because a smaller sound speed enhances the growth of the dark energy fluctuations and this is 
communicated to the matter ones through the gravitational potential. Moreover, the right panel 
shows that for a given sound speed, the non-linear curve is above the linear one because the non- 
linearities enhance the clustering for small scales. This can also be seen in the figures of the previous 
section for the matter perturbations. On the other hand, for k = 0.01/iMpc" 1 (left panel) the blue 
and magenta curves cannot be distinguished because the effect of the sound speed is very small at 
large scales. 

We already commented in the Introduction that the change in the growth function of matter 
induced by dark energy perturbations at the linear level can be at most of order 1% under variations 
of the sound speed. We see from Figure[5]that this is still the case at the non-linear level. However, 
the non-linear corrections on their own have an even somewhat larger effect (~ 3.5%) for small 
scales (0.1 /iMpc" 1 ) and a value of w = —0.8 . 

Notice that we extract the non-linear 5 m from the non-linear matter power spectrum computed 
with the TRG; just by taking the square root of it. This is an approximation with a very small 
error that peaks for small scales and that we can neglect. 




(7.4) 




(7.5) 




(7.6) 
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a a 

Figure 5: Cold dark matter growth g m as a function of the scale factor. It is computed for w = —0.8 at 
k = 0.01 /iMpc -1 (left) and k = 0.1 /iMpc -1 (right) . The black dashed curve represents the smooth linear 
case and the black continuous one the smooth non-linear one. The magenta (c s 2 = 0.001) and the blue 
(0) show the effect of varying the sound speed of dark energy. The continuous blue is non-linear and the 
dashed blue is linear. The dash-dotted magenta is non-linear and the dotted is linear. All the non-linear 
calculations have been done with the TRG. 




0.2 0.4 0.6 0.8 1.0 0.2 0.4 0.6 0.8 1.0 

a a 



Figure 6: Same as Figure \5\ for the total growth g of cold dark matter and dark energy. The green lines 
represent here the ACDM case, that we have included for comparison. The continuous one is the non-linear 
calculation and the dashed one corresponds to the linear result. 



24 



0.45 



0.4C 



w = -0.8 



0.4 



0.6 

a 



0.8 



1.0 



0.45 



0.4C 



w = -0.8 



0.6 

a 



0.8 



1.0 



Figure 7: Cold dark matter growth index 7 m as a function of the scale factor. It is computed for w = —0.8 
at k = 0.01 /iMpc" 1 (left) and k = 0.1 /iMpc" 1 (right) . The black dashed curve represents the smooth linear 
case and the black continuous one the smooth non-linear one. The magenta (c s 2 = 0.001) and the blue 
(0) show the effect of varying the sound speed of dark energy. The continuous blue is non-linear and the 
dashed blue is linear. The dash-dotted magenta is non-linear and the dotted is linear. All the non-linear 
calculations have been done with the TRG. 



The Figure [6] displays the total growth function for cold dark matter and dark energy. The effects 
that we have described for matter also hold in this case. The linear total growth function can be 
easily computed by solving the linear Euler and continuity equations to get the total perturbation 
5tot = Qm<>rn + ^x^x ■ For the non-linear case, let us remind that in our approximation, we actually 
treat 5 X linearly. 

The Figure [6] tells us that the contribution of dark energy perturbations to the total growth is 
relevant. It induces a difference with respect to the ACDM case (where dark energy perturbations 
are absent by definition) that can be as large as ~ 15% for a ~ 0.4 — 0.5 . This result is roughly 
independent on whether the comparison includes the non-linearities or not and the scale at which 
we look. Interestingly, at zero redshift, the effect of changing the sound speed from zero to one is 
also of roughly ~ 15% for k = 0.1 /iMpc" 1 . 

For completeness, we also provide the matter growth index computed from (|7.6p in Figure 0. 
We observe that the effect of the non-linearities is negligible for large (linear) scales and becomes 
of the order of approximately 10 — 15% at small (non-linear) scales. As we can see, this variation 
depends mildly on the redshift (for low values of it: z ~ — 1) . This is an important result 
because it means that non-linearities should be taken into account when using the growth index for 
model discrimination or parameter estimation. Finally, let us remark that this result is essentially 
independent on the sound speed of dark energy, as it can be checked comparing the relative heights 
of the curves in Figure [7]. 
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8 Conclusions 



The near future will likely bring to us a wealth of information, and constraints, on the physical origin 
of the acceleration of the universe, which may be due to dark energy or some modification of general 
relativity. Since at the moment we have almost no clue, neither from observations nor from theory, 
about the possible nature of dark energy, it is a sensible attitude to parametrize it in a general way. 
If dark energy is described as a fluid with a certain equation of state different from — 1 , we can 
study it at the perturbation level using an arbitrary sound speed and anisotropic stress. In this 
work, we focus on the effect of the speed of sound (which must be non-adiabatic) . An immediate 
consequence of a sound speed of dark energy different from the one of light is the emergence of a 
new length scale, the dark energy sound horizon, leading to an effectively scale dependent growth 
function for dark matter, a feature commonly associated to modified gravity scenarios. 

We have explored the consequences of a non-relativistic speed of sound for dark energy, both 
at the linear and the non-linear levels, emphasizing the need of using the rest frame sound speed. 
An accurate treatment of non-linearities is nowadays widely acknowledged to be crucial in order to 
compare theory with future galaxy surveys or measurements of cosmic shear. The traditional third 
order Eulerian perturbation theory is known to fail at redshifts smaller than unity, so that the use of 
resummation procedures, such as the TRG, that can help to extend the validity of semi-analytical 
methods to smaller scales and lower redshifts is mandatory. 

N-body simulations represent a solidly established alternative. However, if one is ultimately 
interested in (for example) parameter forecasting, numerically faster tools, able to scan more effi- 
ciently over a multidimensional parameter space, or over different models, are definitely useful. For 
this kind of analysis, semi-analytical methods (which allow a deeper understanding of non-linear 
effects) are needed. Moreover, in scenarios radically different from ACDM , such as modified gravity 
or massive neutrinos, N-body simulations are still in their infancy. In particular, in the case of 
clustering dark energy with arbitrary sound speed, as far as we know, no N-body simulations have 
been performed yet. Therefore, this work represents the first investigation of non-linearities in these 
scenarios. 

Our main finding is that the dependence of the non-linear correction to the power spectrum on 
the speed of sound is below the percent level in the baryon acoustic oscillation range of scales. This 
is good news for an efficient modelling of the power spectrum in applications such as parameter 
forecasting via, e.g. Markov chains. Indeed, once the effect of the sound speed is taken into account 
at the linear level, the non-linear part can be computed with good accuracy once and for all with 
the model of smooth dark energy with the same equation of state. We have also found that the 
effect of the non-linearities on the matter growth index must be taken into account for small scales 
for a reliable use of this parameter in future applications. 

The observational perspectives for these models are quite challenging but promising. The clean- 
est way to probe this type of scenarios is, in our opinion, through tomographic measurements of 
the gravitational potential via cosmic shear. Indeed, the discriminating power of galaxy surveys is 
probably limited, since the bias issue is further complicated by the emergence of a new clusterized 
component at late times. On the other hand, since lensing measurements typically probe the power 
spectrum at wavenumbers k > 1 — 10/iMpc^ 1 , an accurate assessment of the non-linear effects is 
needed. This paper represents a first step in this direction. 



26 



A A single fluid treatment for zero sound speed 



For the strict case of zero sound speed of dark energy, a treatment of the perturbations in terms of 
just two fields is proposed in [33] . The corresponding equations are 



*r(k) + C(r)0(k) + ydWqMk-p-q)a(q,p)#(q)Mp) =0 (A.l) 
6(k)+W(k) + ^n m H 2 6 T (k) + y>p<i 3 qMk-p-q)/3(q,p)#(q)#(p) =0, (A.2) 



where 5t is related to the total clustering density perturbation 8tot 

5 T = ^ = 6 m + 5 x ^L = 6 -^ (A.3) 

The variable 8 is the common velocity perturbation of matter and dark energy (the same for both 
fluids) and 

C(r) = l + (l + «;)ik. (A.4) 

^ '•m 

It is instructive to compare these equations with our non-linear formalism. Looking at the equation 
(12. 2|) it is not obvious that for a zero sound speed of dark energy (12.221) the velocity divergences of 
dark matter and dark energy are equal. However, this can be seen using the following expression 
for the rest frame sound speed in the Newtonian approximation: 

VP + uP = c 2 (Vp + up) , (A.5) 

which derives from the more general covariant definition [13] : 

(g^ + u»u u ) d v P = c 2 (g^ + tAt") d u p (A.6) 

The equation (|A.5P can be obtained neglecting metric perturbations and combining the time and 
spatial components of (IA.6j) . It can be easily checked that (|A.5|) gives consistently the transforma- 
tion (|2.27p when expanded at first order. Clearly, if c 2 in ()A.5P is zero, the equation (|2.2p reduces 
to that of matter and therefore 8 = 8 X = 8 m , provided that the initial conditions are the same for 
both quantities. 

In reference |43| the dark energy equation for zero sound speed is written as 



6(k) - 3wU8(k) + (1 + w)9(k) + J d 3 p d 3 q5 D (k - p - q)a(q, p)0(q)6(p) = , (A.7) 

which combined with its analogous for w = gives exactly (jA.ip . However, if we set c 2 = in 
(I3.5P we see that the result differs from (|A.7p in a linear 0-term and some 0(2) terms which are 
not included in (|A.7p . We will now explain these differences. 
The continuity equation (12.ip can be rewritten as 



+ 2 u • uP + 3H(p + P) + V • (pu) + P V • u + u • ( VP + uP) = , (A. 



whose last term on the left hand side vanishes for zero rest frame sound speed. In the derivation 
of (jA.ip given in [33] , the terms 2u • iiP and PV • u are both neglected, but it is clear that 
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(at least) the second one of them cannot be discarded because it is of the same order as V • (pu) . 
Actually, making the correct sound speed assignment at first order with (|2.27p and keeping the term 
PV • u , but neglecting 2u • iiP (which gives subdominant contributions) , we obtain the equation 
(]A.7p plus a linear correction in 9 and a non-linear one in 9 x 9 which are both suppressed by 
O {T~L 2 /k 2 ^ . On the other hand, neglecting the term (2u • iiP) in (IA.8j) , amounts at second order 
to take — 2Pu (Hu + V0jv) = , as it can be directly read from (|2.5p or (|2.2p . This approximation 
means setting to zero several second order non-linear contributions to the perturbed continuity 
equation that are proportional to w (and hence vanish for dust but not for dark energy). With 
some algebra, one can check that these corrections are also suppressed at small scales by at least 
one power of 7i 2 /k 2 (where k refers to any non-linear combination of momenta) with respect to 
the usual non-linear a term . In conclusion, in order to get (lA.7h and hence (lA.ip we just have to 
discard linear and non-linear terms in the perturbed continuity equation for dark energy that are 
small at scales k S> H ■ Notice, incidentally, that to compute the precise form of these negligible 
second order corrections we would need to know (I2.27|) at second order; but, however, the linear 
one coming from PV • u can be readily obtained (and it is a term —9w(l + w)H 2 9/k 2 on the right 
hand side of (|A.7j) ) . 

Notice that if we take (|6.2p and replace and 5^ by their non-linear counterparts, we obtain 
precisely (|A.2j) . The equations (jA.ip and (IA.2P are a good approximation for the non-linear evo- 
lution of the total density perturbation in the case of dark energy with zero sound speed, as the 
results of Section [7] show. This approximation is used in [59] to apply the TRG for the computation 
of the total power spectrum in the case c 2 = of maximal clustering Q. 
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